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[i] This study adopts a “precipitation object” approach by using 14 years of Tropical 
Rainfall Measuring Mission (TRMM) Precipitation Feature (PF) and National Centers for 
Environmental Prediction (NCEP) reanalysis data to study rainfall structure and 
environmental factors associated with extreme heavy rain events. Characteristics of 
instantaneous extreme volumetric PFs are examined and compared to those of intermediate 
and small systems. It is found that instantaneous PFs exhibit a much wider scale range 
compared to the daily gridded precipitation accumulation range. The top 1% of the rainiest 
PFs contribute over 55% of total rainfall and have 2 orders of rain volume magnitude greater 
than those of the median PFs. We find a threshold near the top 10% beyond which the PFs 
grow exponentially into larger, deeper, and colder rain systems. NCEP reanalyses show that 
midlevel relative humidity and total precipitable water increase steadily with increasingly 
larger PFs, along with a rapid increase of 500 hPa upward vertical velocity beyond the top 
1 0%. This provides the necessary moisture convergence to amplify and sustain the extreme 
events. The rapid increase in vertical motion is associated with the release of convective 
available potential energy (CAPE) in mature systems, as is evident in the increase in CAPE 
of PFs up to 10% and the subsequent dropoff. The study illustrates distinct stages in the 
development of an extreme rainfall event including (1) a systematic buildup in large-scale 
temperature and moisture, (2) a rapid change in rain structure, (3) explosive growth of the 
PF size, and (4) a release of CAPE before the demise of the event. 

Citation: Zhou, Y., W. K. M. Lau, and C. Liu (2013), Rain characteristics and large-scale environments of precipitation 
objects with extreme rain volumes from TRMM observations, J. Geophvs. Res. Atmos., 118, 9673-9689, doi: 10. 1002/jgrd. 50776. 


1. Introduction 

[ 2 ] One of the most likely consequences of global wanning 
is an increase in extreme weather events, such as droughts, 
floods, and heat waves [IPCC, 2012]. A number of recent 
studies have shown increased frequency and intensity of 
extreme precipitation events in recent decades [ Easterling 
et al., 2000; Trenberth et ah, 2003; Groismcm et al., 2005; 
Lau and Wu, 2007; Allan and Soden, 2008; Min et al., 
2011]. Unlike the global mean precipitation, which is 
constrained by the global radiative convective balance and 
increases at a rate (~3%°C _1 ) much slower than the atmo- 
spheric water vapor (~7%°C _1 ) [Allen and Ingram, 2002; 
Held and Soden, 2006], extreme precipitation is more 
closely related to moisture convergence in the atmosphere 
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and increases at a much faster rate than the mean precipita- 
tion [Kharin et al., 2007; O’Gorman and Schneider, 2009; 
Lau and Wu, 201 1; Lau et al., 2013]. 

[ 3 ] However, it is a challenge to pin down the direct con- 
nection of an individual extreme precipitation event to global 
wanning. Extreme precipitation occurs under a variety of 
mesoscale and synoptic meteorological conditions. Many 
extreme precipitation events can be identified with mesoscale 
convective systems (MCSs), such as tropical storms and hur- 
ricanes [Shepherd et al., 2007; Lau et al., 2008], squall lines, 
the organization of mesoscale convective complex (MCC), 
and heavy-precipitation supercells that often produce flash 
floods [Moore et al., 2003; Smith et al., 2001]. Extreme pre- 
cipitation can also result from terrain-forced convection 
[Petersen et al., 1999; Pontrelli et al., 1999] and from both 
strongly and weakly forced synoptic systems [Maddox 
et al., 1979; Heideman and Fritsch, 1988; Bradley and 
Smith, 1994], Despite this, extreme precipitation is usually 
a result of many coincident meteorological conditions, such 
as a buildup of large local atmospheric instability or persis- 
tent large-scale moisture transport and convergence [Reale 
et al., 2012], Orographic forcing plays a significant role in 
many cases, such as the disastrous flooding caused by 
Typhoon Morakot in 2009 in Taiwan [Ge et al., 2010; Lee 
et al., 2011] and rainstorms in 2010 in Pakistan [Houze 
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et al., 201 1], The interaction of the storm's large-scale envi- 
ronment and life cycle is very important and unique to each 
extreme event and deserves careful study. However, detailed 
study is impossible for every extreme event because of the 
lack of thermodynamic observations of the entire stonn life 
cycle and the large number of such events. Moreover, even 
a detailed case study may not be able to illuminate the 
connection of individual events to changes in large-scale 
environments that are attributable to global warming, as in 
the case of the Pakistan flood in 2010 [Wang et al., 2011; 
Lau and Kim, 2012; Dole et al., 201 1] and Typhoon Morakot 
[Lee et al., 2011], 

[ 4 ] Therefore, regardless of whether or not the cause is 
global warming, extreme events are always related to changes 
in the large-scale environment and local feedback processes 
that amplify their development. In this study, we use a 
composite technique to examine the basic characteristics of 
extreme precipitation and the associated mean large-scale en- 
vironment in order to gain insight on the changing rain struc- 
ture and the environmental conditions that are more likely to 
produce extreme events. Section 2 describes the data set and 
methodology. Section 3 shows rain characteristics for different 
sizes of precipitation objects. Section 4 examines mean large- 
scale environments associated with extreme precipitation 
events. A summary and discussion is provided in section 5. 

2. Data and Methodology 

[ 5 ] The Tropical Rainfall Measuring Mission (TRMM) 
has accumulated more than 14 years of state-of-the-art tropi- 
cal rainfall products including rain rates, types, structures, 
latent heating profiles, and related water vapor and cloud 
properties. These data products offer a unique opportunity 
for research on extreme precipitation characteristics and their 
modulation by climate processes. Here we use the University 
of Utah Precipitation Feature (PF) database [Liu et al., 2008], 
which provides rainfall characteristics for individual precipi- 
tation objects defined by grouping contiguous satellite rainy 
pixels with related cloud and precipitation observations. 
This database was developed within the framework of PFs 
earlier defined by Nesbitt et al. [2000] and was updated with 
version 7 of TRMM products. The characteristics of each PF 
event, which include the raining area, total rain volume, 
cloud top temperature inferred from infrared brightness 
temperature, precipitation radar echo top height, lightning 
flash rates, etc., are derived from collocated measurements 
and retrievals from the Precipitation Radar (PR), TRMM 
Microwave Imager (TMI), Visible and Infrared Sensor, and 
Lightning Imaging System. This database covers the entire 
TRMM period (1998 to present) from 36°S to 36°N, enabling 
quick identification of specific cloud and precipitation sys- 
tems and characterization of these systems from multiple 
TRMM instruments. The surface meteorology and vertical 
profiles of temperature, geopotential height, wind, and humid- 
ity associated with each PF were extracted from associated 
six hourly 2.5° * 2.5° National Centers for Environmental 
Prediction (NCEP) reanalysis [Liu et al, 2008; Kalnay et al., 
1996; Kistler et al., 2001]. The NCEP data were first tempo- 
rally interpolated to the PF time, and the nearest grid point to 
the center of the PF was selected for each PF. In this study 
we use a level-2 radar-based PF data set that was defined by 
grouping contiguous raining pixels indicated by PR retrievals 


[Iguchi et al., 2000, 2009], Each PF in the data set has at 
least four continuous PR pixels with nonzero surface rain- 
fall to remove detection noise. We focus on the PFs with 
extreme values of volumetric rainfall from PR retrievals 
and associate them with the large-scale environments from 
collocated NCEP reanalysis data. 

[6] The precipitation feature database is similar to the cloud 
object approach developed at NASA Langley Research Center 
to study the cloud and radiation characteristics of various 
cloud objects, i.e., deep convective clouds or shallow clouds 
[Xu et al., 2005, 2007; Eitzen and Xu, 2005]. Composite 
and clustering techniques were used in the cloud object 
studies to categorize cloud properties in different types of 
clouds. Lin et al. [2006, 2007] have studied the effect of 
environmental conditions on tropical deep convective sys- 
tems observed from the TRMM satellite and have found that 
not only rainfall rate but also rainfall efficiency increases 
with sea surface temperature. Although different names 
are used, both groups use an object-oriented approach that 
treats continuous satellite footprints of cloud or precipita- 
tion areas as a single object with physical and statistical 
properties. Our study based on PF data could also be viewed 
as a study of precipitation objects. 

[ 7 ] In addition to the PF data, we use the TRMM 3B42 
[Huffman et al., 2007] daily rainfall product to provide addi- 
tional information for the time-accumulated rainfall informa- 
tion for the grid box where PF locates. The 3B42 algorithm 
merges precipitation estimates from passive microwave 
(PMW) products from TRMM TMI and Special Sensor 
Microwave Imager (SSM/I) and rain estimates from geosta- 
tionary infrared (IR) products. The SSM/I data are calibrated 
to the TMI measurements, while the IR estimates are cali- 
brated with PMW estimates. The daily 3B42 data are aggre- 
gated from the original three-hourly temporal resolution with 
0.25° x 0.25° spatial resolution covering a latitudinal domain 
from 50°S to 50°N. 

[s] The definition of extreme precipitation is highly depen- 
dent on the application and data set available. Depending on 
applications, various indices have been used, e.g., maximum 
1 or 5 day precipitation, daily intensity index, or the number 
of heavy and consecutive precipitation days [Zhang et al., 
2011]. Different authors define extreme thresholds based 
on different percentiles of the precipitation distributions, 
with the farther tail indicating more extreme events. In this 
study, we define extreme events using the volumetric rain 
rate (VRR: a product of instantaneous rain intensity and 
the size of the rain system) from PR near-surface rain 
[Iguchi et al., 2000, 2009] of PFs to focus on the total rain 
amount the system produces at the TRMM overpass time. 
VRR pertains to an instantaneous and intrinsic measure of 
the rainfall intensity that approximates total latent heat 
release of the system at any instant. Using VRR as a mea- 
sure of rain intensity, we gain an in-depth perspective of 
the fundamental physical processes governing the occur- 
rence of extreme precipitation events. The correlation of 
extreme PF with daily extreme rainfall from TRMM 3B42 
data will be discussed in section 3.2. 

[ 9 ] We also point out that extreme PFs defined based on 
VRR are focused on stonn extent in size and rain amount 
rather than stonn convective intensity, which is characterized 
by strong updrafts, high lightning rates, and elevated radar 
reflectivity [Zipser et al., 2006], 
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Rain percentile 



Rain percentile 

Figure 1. (a) Volumetric rain thresholds (left y axis, line with diamond) and percentage of cumulative rain 
amount (right y axis, line with cross) with percentiles in rain volume. The percentile represents the percent- 
age of PFs above the threshold over all PFs in 36°S-36°N during 1998-201 1 . (b) Similar plot as Figure la 
but using nonzero 0.25° x 0.25° grid 3B42 data in 36°S-36°N during 1998-2011. 


3. Results 

3.1. Probability Distributions of PF Rainfalls 

[to] One advantage of using PF VRR as a reference for 
extreme rain events is the large data sample based on satellite 
pixels. There are more than 28 million PFs in the database 
from January 1998 to December 2011 with a minimum of 4 
PR pixels or 75 km 2 in size. The smaller systems (less than 
4 PR pixels) excluded from the PF database are not expected 
to contribute to a significant portion of the total rainfall, much 
less the top 1% of the extreme events [Liu, 201 1]. To begin, 
we computed the probability distribution function (PDF) of 
all PFs according to their VRR and determined the thresholds 
for every 10% in frequency, as well as the top 1%, 0.1%, 
0.01%, and 0.001% for extreme events (Figure la). The 
thresholds increase nonlinearly as the percentile tends toward 
extreme categories, with an abrupt jump near the top 1 0th 
percentile. For example, VRR is less than 2300mm/hkm 2 
below the top 10th percentile but reaches above 40,000 and 
575,000 mm/h km 2 in the top 1% and 0.01%, respectively, 
which is more than 2-3 orders greater than that of median 
PFs. On the other hand, Figure la shows that PFs from the 
top 1% contribute to over 55% of the total instantaneous rain- 
fall, while those below the 10th percentile (or the lower 90%) 
contribute less than 20% of the total rainfall. It should be 
noted that even though the majority (70-80%) of PFs in 
Figure la are from oceans, the PDFs and the corresponding 
thresholds from land and oceans separately are extremely 
close (figure not shown), so only one set of thresholds is used 
in this study. 

[t l] Since both rain rate and volumetric rain rate of PF are 
snapshot measurements, it is interesting to know how those 
snapshots translate into the time-integrated rain rate. The daily 
3B42 data from TRMM [Huffman et al., 2007] contain 
enough samples (some microwave overpasses and continuous 


IR) to represent daily rain accumulation fairly well. Physically 
speaking, a large rainy system usually lives for a long time and 
brings a heavy rainfall accumulation locally. Therefore, a 
snapshot of large rainy systems from TRMM would be very 
likely to lead to a heavy accumulation event in daily precipita- 
tion. Flowever, a small system with an intense rain rate locally 
does not necessarily lead to a heavy rainfall accumulation due 
to a potentially short lifetime. On the other hand, a heavy daily 
accumulation event can be from a relatively small system last- 
ing a long time period due to continuous moisture supply, e.g., 
orographic lifting of monsoonal flow. In order to provide con- 
text and contrast PF VRR data with time-integrated precipita- 
tion data, we show similar statistical and spatial distribution 
plots for the daily 3B42 data. 

[ 12 ] Figure lb shows a more linear curve of the PDF forthe 
daily 3B42 data than the PF data. The top 1% of the daily 
gridded precipitation contributes only 13% of the total rainfall 
versus 55% in the case of the instantaneous PF. Compared to 
the 3B42 time-averaged gridded data, the PF-based PDF 
shows that instantaneous rainfall has a larger scale range, with 
a much longer tail at the extreme end of the distribution. 
Despite the low frequencies for those above the top 1%, the 
absolute number of extreme PFs is still quite large due to the 
enormous total PFs globally (Figure 2), providing a sound 
basis for our statistical analysis. 

3.2. Geographical Distributions of Extreme PFs 

[ 13 ] The geographical distributions of the most extreme 
rainy PFs in VRR (Figure 2, top) differ significantly from 
those of the most intense storms [Zipser et al., 2006, 
Figure 3] and most intense daily rainfall (Figure 2, bottom). 
The most intense storm was defined by the strength of con- 
vection; therefore, elevated high reflectivity, high lightning 
rate, and strong ice signature at microwave channels (i.e., 
minimum 85PCT) were used as the main criteria to identify 
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Figure 2. (top) Geographical distribution of PFs with extreme volumetric rain. Contours are for popula- 
tion of top 1% of extreme PFs in 2° x 2° grid boxes. TRMM sampling bias has been removed in the 
population contours by dividing the ratio of actual and mean sample numbers in the grid box. Green, 
purple, and black triangles are the locations of the top 0.1%, 0.01%, and 0.001% extreme PFs. Note that 
there are more triangles at subtropics partially due to the sample biases near 30° latitudes, (bottom) 
Number of top 1% events in 2° x 2° grid boxes (> 72 mm/d) from daily 3B42 data during 1998-2011. 


the most intense storms. Figure 2 (top) shows that most of 
the extreme rainy PFs in VRR occur over the oceans as 
opposed to over the land regions in the case of severe 
storms. High concentrations of systems with extreme volu- 
metric rain can be seen over the Northwest Pacific, South 
Pacific Convergence Zone (SPCZ), Northwest Atlantic, 
and the Bay of Bengal, accompanied by moderate concen- 
trations along the tropical Intertropical Convergence Zone 
(ITCZ). Over land, large volumetric PFs are found along 
the coastal South China Sea and the southeast U.S. as a 
result of tropical storms moving in from the oceans and 
large moisture transport from the sea. High concentrations 
of extreme PFs are also found over South America south 
of the Tropic of Capricorn — the same region frequented 
by intense stonns [Zipser et al., 2006]. 

[ 14 ] It is important to understand the relationship between 
extreme PFs in VRR and extreme precipitation accumulation 
events at the surface, since the former only captures an in- 
stantaneous snapshot while the latter requires rainfall over 
a certain duration. Figure 2 (bottom) shows the geographi- 
cal distribution of the top 1% of extreme daily events (daily 
precipitation greater than 72 mm/d) based on the 3B42 data. 
The most frequent occurrences of daily extreme rainfall are 
found over the northwest Pacific and South China Sea 
along the tropical typhoon route, SPCZ, the Pacific and 
eastern Atlantic ITCZ, the Bay of Bengal, Gulf of Mexico, 
Southwest U.S., Central Amazon, northeast Argentina, 
and the Gulf of Guinea. It is obvious that regions with high 
concentrations of extreme PFs also have a high frequency 
of daily extreme events. However, Figure 2 (top) shows 
extreme PF events quite evenly spread out in the equatorial 
Indian Ocean, west Pacific, and SPCZ region (even for the 
very extreme ones), while the 3B42 map (Figure 2, bottom) 


shows regions of high concentration of extreme events at 
the center of SPCZ and equatorial ITCZ. Also from 
Table 1, many extreme PFs do not produce extreme daily 
rainfall. This evidence shows that many extreme PFs do 
not last long enough to produce large rain accumulation. 
The meteorological conditions used to produce extreme 
PF quickly dissolve. This indicates that an extreme daily 
precipitation accumulation event might require a more 
rigorous environment and higher level of organization of 
the precipitation system. Particularly worth noticing is the 
high frequency of daily extreme rainfall off the west coast 
of India and the southeast coast of Asia as a result of the 
orographic effect. 

[ 15 ] To obtain a more quantitative measure, we collocated 
top PFs and their corresponding daily mean precipitations 
from 3B42 and calculated the percentage of top PFs that pro- 
duced daily extreme rainfall (Table 1 ). It was found that only 
about 17% of PFs in the top 1% were associated with heavy 
daily events (greater than 50 mm/d). This rises to 36.4% and 
54.9% for the top 0.1% and 0.01% of PFs, respectively. Since 
many of the large rain systems span over night, we also com- 
puted 3 day rain accumulation centered with the PFs. The 
percentage of cases with 3 day rain accumulation greater than 
50 mm rises to 40.3%, 61.5%, and 77.2% for PFs in the top 
1%, 0.1%, and 0.01%, respectively, which is consistently 
larger than that of 1 day accumulation. These results indicate 
that larger PFs usually last longer, as shown in Liu [2011], 
and many of the extreme PFs extend over more than 1 day. 
The above results show that extreme PFs (especially for those 
above top 0.1%) are most likely associated with large accu- 
mulated rain events, so the results shown in this study will 
generally apply to extreme events in time-averaged measure- 
ments (a day or shorter). Of course, there are instances of 
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Figure 3. Normalized total number of top 1% extreme PFs (shading) in 2° x 2° grid boxes and mean pre- 
cipitation from 3B42 data (contour) in each season. TRMM sampling bias has been removed by dividing 
the ratio of actual and mean sample numbers in the grid box. 
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large instantaneous rain showers from short-lived systems, 
such as afternoon showers in deep tropical land areas, that 
are captured by extreme PFs but do not produce heavy daily 
rainfall [Takavabu et al., 2012], On the other hand, heavy 
daily rainfalls could be the result of persistent moderate rain 
throughout the day [Liu, 2011], 

3.3. MCS and Extreme PF 

[i6] While it is likely that extreme volumetric PFs are also 
large in size (see Figure 5 for more detailed discussion), here 
we will examine what kind of rain systems (e.g., MCS, warm 
rain system, or others) are found in different percentiles of 
PF. MCS is one of the dominant heavy-precipitation systems 
in the tropics, with complex organized storms that extend 
more than 100 km in one direction but are smaller than 
extratropical cyclones [Houze and Churchill, 1987]. Tropical 
cyclones are one of such systems that contribute more than 
60% of the total extreme events in northwest Pacific based 
on pentad Global Precipitation Climatology Program data 
[Lau et al., 2008]. Other MCSs include MCC, squall lines, 


lake-effect snow events, and polar lows. Based on 1 year of 
SSM/I data, Mohr et al. [1999] found that MCSs consisted 
about 10%-20% of the sampled systems by number but con- 
tributed 70%-80% of the rainfall throughout the Tropics. 
Using a similar MCS definition (PF size > 2000 km 2 ) as 
Mohr et al. [1999], we find that very few PFs below the top 
10th percentile are contributed by MCS simply due to the size 
requirement by the definition (Table 2). About 79.6% and 
67.8% of the PFs in the top 1% have one or several MCSs 
imbedded in the system over the land and the oceans, respec- 
tively. For those at the top 0.1%, a larger percentage (90%) 


Table 1. Percentage of Top PFs With 1 and 3 Day Precipitation 
Accumulation Centered at PF Locations and Times From TRMM 
3B42 Greater Than 50 mm/d 



1 Day Accumulation 

3 Day Accumulation 


>50 mm 

>50mm 

Top 1% 

17.0% 

40.3% 

Top 0.1% 

36.4% 

61.5% 

Top 0.01% 

54.9% 

77.2% 
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Table 2. Fraction of PFs Categorized by Volumetric Rainfall That 
Are Identified as MCSs a 



<10% 

10% 

1% 

0.1% 

0.01% 

0.001% 

Land 

0.0 

5.2 

82.7 

93.9 

98.4 

100.0 

Ocean 

0.0 

2.7 

72.0 

91.1 

96.7 

99.6 

Global 

0.0 

3.5 

74.8 

91.7 

97.0 

99.7 


a With size > 2000 km 2 . 


could be characterized with MCS globally, with an even 
higher percentage over some land regions; i.e., 99% of the 
PFs in the top 0.1% are identified as MCSs in the Bay of 
Bengal and in equatorial West Africa. These results indicate 
a much lower percentage (about 3 ~4%) of MCS in all PFs 
compared with the results of Mohr et al. [1999]. There are 
a few explanations for this discrepancy. The SSM/I single 
pixel size is 13 x 15 km 2 , larger than the smallest PF 
(75 km 2 ) in this study, which could lead to some small 
systems missing from Mohr et al. [1999] that are detectable 
by PR. More importantly, however, studies using SSM/I 
rely on the ice scattering signal, which detects the deep 
convection with ice but totally misses warm rain systems 
as a result. Therefore, the total population of systems in 
Mohr et al. [1999] would be smaller due to missing warm 
rain systems, especially over the oceans. 

[n] Even though it has been shown that most warm rain 
systems are small in size [Liu and Zipser, 2009], it is not clear 
what percentage of the PFs in each volumetric rain category 
belongs to warn rain systems or whether warm rain can also 
produce extreme rainfall. Defining warm rain events as those 
that have a minimum IR temperature greater than 273 K, a 
minimum 85 GHz PCT (polarization-corrected brightness 
temperature) greater than 250 K, and a maximum storm 
height less than 4.5 km, we find that warm rain events con- 
tribute only 1.4% and 0.2% to the top 10th percentile over 
the oceans and the land, respectively, which means that wann 
rain events rarely contribute to extreme PFs by rain volume 
(Table 3). Interestingly, only about one third of the PFs in 
the 90th-30th percentiles belong to warm rain systems, and 
an even smaller percentage of PFs in the bottom 10% meet 
the warm rain definition. Therefore, it is likely that a large 
number of small size and low volumetric PFs are either 
isolated small convective cells or patches of PF separated 
from a large system. On the other hand, a considerable num- 
ber of PFs (~20% in the top 1% and 10% in the top 0.1%) are 
of neither MCS nor pure warm rain systems but most likely 
have mixed features of warm or mixed-phase rain attached 
to a colder system [Liu and Zipser, 2009]. We want to point 
out that the minimum size of 4 PR pixels (75 km 2 ) for this 
study excludes a large number of small warm rain systems, 
because most of the warm rain systems are less than 30 km 2 
in size [Liu and Zipser, 2009], Nevertheless, it does not 
change the conclusion that warm rain contributes very little 
to extreme volumetric rain events. 

3.4. Seasonal and Diurnal Distributions of Extreme PF 

[is] Seasonal variations of extreme precipitation are similar 
to seasonal variations of mean precipitation (Figure 3). For ex- 
ample, high concentrations of extreme events are found in the 
maritime continent and Southeast Asia in June-July-August 
(JJA) and are associated with the monsoon system. Extreme 
precipitation associated with the tropical ITCZ and SPCZ is 


also in proximity to the intensity and spread of the monsoon 
system. There is a much broader region of heavy mean precip- 
itation and frequent occurrences of extreme PF events in the 
SPCZ in December-January-February (DJF) and March- 
April-May than JJA and September-October-November. 
On the other hand, both the East Pacific and Atlantic ITCZ 
regions show more extreme PFs in JJA than in DJF. Over 
land, more extreme events are found over tropical and sub- 
tropical South Africa and the Amazon in DJF than in JJA. 
The seasonal variations of extreme events poleward of 30° 
in both hemispheres are likely associated with seasonal 
migration of midlatitude storm tracks. For example, more 
extreme PFs are found north of 30°N in the Pacific Ocean 
in the winter than in the summer due to an equatorward 
migration of the jet stream. These storms are distinct from 
the tropical systems in the maritime continent and Southeast 
China Sea, where an elevated number of extreme PFs appear 
in JJA due to monsoon and tropical cyclones. 

[ 19 ] To illustrate the diurnal distributions of PFs according 
to their volumetric rainfall, we computed normalized fre- 
quencies for each hour for PFs within 10%, between 10% 
and 0.1%, and larger than 0.1% in four subregions: tropical 
land and oceans (20°S-20°N) and subtropical land and 
oceans (20°N-36°N and 20°S-36°S) as shown in Figure 4. 
A relatively weak diurnal cycle is found over the oceans 
compared to over land, as shown by many previous studies 
[Nesbitt and Zipser, 2003; Bowman et al., 2005; Liu and 
Zipser, 2008; Liu, 2011], Because land has low heat capacity 
compared to the ocean, it has a much larger diurnal cycle of 
surface temperature than the ocean. The diurnal cycle of 
PFs over land is mostly driven by convective overturning 
caused by afternoon heating of the surface and the lower 
boundary layer. However, the oceanic diurnal cycle is con- 
sidered to be related to the much weaker daytime stabiliza- 
tion of solar heating and night destabilization of radiative 
cooling at the cloud top [Kraus, 1963]. For tropical oceans, 
smaller systems (<10% and 1 0%— 0. 1 % categories) show 
weaker diurnal variations than the larger systems (> 0.1%). 
They also reach peak frequency earlier in the morning (3 A.M.) 
than the larger systems (6-9 A.M.). This is consistent with 
Nesbitt and Zipser [2003], who showed that small systems 
peak at 3 A.M. but the number of MCSs peaks at 7 A.M. 
over the oceans. The largest diurnal variation occurs for 
the top 0.1% of PFs; this is probably due to continuous 
expanding of smaller MCS during the early morning 
before dying down in the afternoon [Nesbitt and Zipser, 
2003], Over the tropical land, the largest diurnal variation 
occurs in the top 1 0%— 1 %, in which surface heating favors 
medium size afternoon convection. The afternoon peaks are 
reversed for the extreme PFs (top 0.1%), in which most 
MCSs prolong into midnight and keep growing until the 
early morning. The diurnal cycle in subtropical land and 


Table 3. Fraction of PFs Categorized by Volumetric Rainfall That 
Are Identified as Warm Rain Systems' 1 



<90% 

90-30% 

30% 

20% 

10% 

>1% 

Land 

9.5 

34.1 

16.3 

7.4 

0.2 

0.0 

Ocean 

12.3 

38.6 

19.5 

9.3 

1.5 

0.0 

Global 

2.0 

9.2 

3.9 

1.6 

1.1 

0.0 


“Minimum infrared brightness temperature > 273 K. 
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Figure 4. Diurnal distributions of the total number of PFs from three PF categories for four subregions: 
tropical (20°S-20°N) oceans, tropical land, subtropical (20°N-36°N and 20°S-36°S) oceans, and subtrop- 
ical land. 


oceans is similar to the diurnal cycle in the tropics with a 
relatively smaller magnitude, which means that the differ- 
ence between the diurnal cycles is mainly due to thermal 
conditions of the surface between land and ocean. These 
results are consistent with various previous studies [ Yang 
and Slingo, 2001; Nesbitt and Zipser, 2003; Liu and Zipser, 
2008; Liu, 2011], 

3.5. Variations of Rain Size and Rain Rate With VRR 

[ 20 ] In this subsection, we investigate the relationship 
between rain size, volume, and mean intensity. In general, 
total volumetric rain is expected to increase with PF size 
(Figures 5a and 5b). Flowever, the relationship between 
VRR and rain size is far from linear, as demonstrated by 
the broad spread of the PDF (Figures 5 a and 5b) and the dif- 
ference of PDFs in difference subregions (Figures 5c and 5d). 
For cases below the median VRR, the PF stays predomi- 
nantly smaller than 150 km 2 and is associated with primarily 
meso-y or storm-scale systems. Above the median VRR, the 
size increases rapidly to 1500-5000 km 2 . Beyond the top 
10th percentile, the growth is exponential, with the size 
increasing 20 times from 5000 to 1 00,000 mm/h km 2 and 
VRR increasing 2-3 orders (Figure 1). The exponential 
growth is more due to the rapid increase in size than in the 
rain rate, as is evident in the near-vertical orientation of the 
PDF beyond the top 10th percentile. 

[ 21 ] As seen from Figures 5a and 5b, tropical oceans dis- 
play a much higher frequency of small PFs (under 150 km 2 ) 


and moderate VRR (80th-30th percentile) compared to trop- 
ical land. Tropical land has a slightly higher frequency out- 
side this core distribution area, indicating more spread-out 
and diversified rain systems over the land regions. The higher 
frequency in the core distribution area in tropical oceans, on 
the other hand, indicates more abundant shallow convection 
that is limited in size and rain production. Compared to the 
PDFs in the tropics, the PDFs in subtropics show increased 
and decreased frequency on the left and right sides of the 
core distribution area, respectively (Figures 5c and 5d). 
This indicates that PFs in the subtropics have a slightly 
higher chance to be larger in size in the same rain percen- 
tiles, possibly due to relatively lower moisture content in 
the subtropics than in the tropics. There could be many 
reasons for large systems to occur in the relatively dry sub- 
tropics, i.e., the wave intrusion of cold fronts from higher 
latitudes and the organization of large MCSs under subtrop- 
ical monsoonal flow (e.g., the Mei-yu front in China, MCSs 
over India and Argentina). The warm ocean currents in the 
subtropics may also contribute to the organization of large 
systems over oceans. Some mesoscale systems, such as 
hurricanes, are organized more easily off the equator due 
to the presence of the Coriolis force. 

[ 22 ] The relationship between PF size and VRR is further 
illustrated in terms of mean conditional rain rates, defined 
as the ratio of VRR to PF size (Figure 6a). Surprisingly, mean 
conditional rain rates from all four subregions are rather close 
to each other and increase gradually from less than 5 mm/h in 
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Figure 5. Two-dimensional histogram ofPFs with size and total volumetric rain for (a) tropical ocean and 
(b) tropical land; (c) difference between subtropical ocean and tropical ocean; (d) same as Figure 5c but for 
land. Dashed lines indicate locations of maximum frequency. 


small PFs (< top 10th percentile) to around 7mm/h in ex- 
treme PFs. However, mean rain rates in convective and strat- 
iform areas have quite different characteristics. The mean 
rain rates in convective areas increase from around 1 mm/h 
in PFs with low VRR to around 4-6 mm/h in the top 20th 
percentile and then increase rapidly to nearly 20 mm/h in 
the largest PFs (Figure 6a, blue lines). Compared to other 
regions, tropical oceans have relatively lower mean convec- 
tive rain rates for PFs up to 1% but relatively higher mean 
convective rain rate for PFs larger than 1%. The mean rain 
rates in the stratiform areas are much smaller than those of 
the convective areas, starting from around 0.6 mm/h in the 
smallest PFs to around 5 mm/h in the largest PFs (Figure 6a, 
red lines). Previous studies have found that the fraction of 
convective rain area is a major factor in determining the mean 
conditional rain rate [. Nesbitt and Zipser, 2003], Figure 6b 
shows larger differences in convective rain area for smaller 
PFs ranging from 20% in subtropical land to 60% in tropical 
oceans. The differences become smaller as VRR increases, 
and convective areas in all regions except subtropical land 
drop to less than 15% as small systems develop into large 
organized systems with extended stratiform rain area. As a 
result, mean conditional rain rates are dominated by stratiform 
rain area for the large PFs, while they are not far apart from 
both convective and stratiform rain rates for small PFs. Even 
though there is a rapid increase of convective rain rate and a 
rapid decrease of convective rain area at 1 0%— 1 %, a much 
smoother change of the mean conditional rain rate across 
different PF sizes is found (Figure 6a, black lines). 


[ 23 ] Besides the reverse trend in convective rain rate and 
area with VRR for all regions, there are regional compensat- 
ing effects which further reduce the differences in conditional 
rain rate. For example, tropical oceans have the smallest con- 
vective rain rates and largest convective rain areas for smaller 
PFs but the largest convective rain rates and next to smallest 
convective rain areas for larger PFs. Subtropical land is the 
opposite of tropical oceans, having the smallest convective 
rain areas and the largest convective rain rates for smaller 
PFs but the largest convective rain areas and the smallest 
convective rain rates for larger PFs. Our results show that 
even with the decrease of fractional convective rain area, 
the conditional rain rates for both convective areas and strat- 
iform areas increase with the PF size, which contribute to an 
overall increase in the mean total conditional rain rates across 
different PF sizes (Figure 6a). This is contradictory to Nesbitt 
and Zipser [2003], who showed that mean conditional rain 
rates in convective areas and stratiform areas are nearly 
constant throughout the day, while the change in convective 
fraction contributes to the diurnal cycle of mean conditional 
rain rate in different rain systems. 

3.6. Variation of Rain Structure With VRR 

[ 24 ] In this subsection, we further investigate the evolution 
of rain characteristics as a function of rain volume. To illus- 
trate the significant departure of extreme PFs from the lower 
99% of the rain systems, we show the changes of mean 
parameters that describe the vertical extension and rain inten- 
sity of the PF with VRR (Figure 7). Figure 7a shows that the 
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Figure 6. (a) Mean conditional rain rate for all (black), convective (blue), and stratiform (red) rain areas. 
The y axis is in logarithmic scale for PFs smaller than 10% and in linear scale for PFs larger than 10%. (b) 
Percentages of convective rain area. Line type indicates four different subregions. 


mean minimum IR temperature (MinlR) separates into two 
groups (land and ocean) for the small PFs at around 250 
and 265 K, respectively, with no latitudinal delineation. 
This indicates that smaller PFs over the ocean reach lower al- 
titude compared to those over land (see also Figures 8a and 
8b). The MinlR drops to below 220 K at the 10th percentile, 
and the difference between the tropics and subtropics starts to 
emerge in place of land and ocean, with MinlR in the tropics 
being 10-20K colder than MinlR in the subtropics due to 
much higher and colder tropopause in the former. The mean 
maximum height (Maxht, Figure 7b) and mean maximum 
40dBZ (Maxht40, Figure 7c) also show a dramatic increase 
of storm height as VRR reaches the 10% — 1% range. The 
mean storm height for small PFs also delineates by land 
and ocean at 3.5 and 5.5 km, respectively, but tend to merge 
as all of them rise to above 15 km for the most extreme PFs. 
Most small PFs are too weak to produce a 40dBZ radar 
signal, so the mean Maxht40 is near zero until the 40th 
percentile but rises rather dramatically at 1 0% — 1 % to above 
4 km for the most extreme PFs. 

[ 25 ] The 85 GHz PCT reflects scattering signals of ice 
hydrometeors (graupel and hail) from convective updrafts 
above the freezing level [Spencer et al., 1989], The minimum 
37 GHz PCT is sensitive to larger hydrometers and can be 
used to identify low-level rain systems that are obscured by 
the precipitation-sized ice particles aloft. Both minimum 
85 and 37 GHz PCTs show significant depression for the 
extreme PFs starting at the 10th percentile, indicating deep 
convection and strong ice scattering signals from the larger 
systems (Figures 7b and 7d). Lightning flash rates are highly 
correlated with convective ice water paths [Petersen and 
Rutledge, 1998; Petersen et al., 2005], The mean flash rate 
also begins to increase at the 10th percentile, with much 
higher flash rates over the land than over the oceans for the 


PFs with similar VRR (Figure 7c). In all the cases, PFs over 
land have stronger convective intensity proxies than those 
over the oceans for the same VRR, consistent with much 
stronger convection over land. All parameters examined show 
dramatic change when VRR reaches a threshold of the 10th 
percentile (2300mm/hkm 2 ). 

[ 26 ] Even though the mean characteristics of PFs show a 
systematic change with VRR, a large standard deviation 
is expected for any parameter at any given rain range. 
Figure 8 shows a two-dimensional histogram of Maxht 
versus MinlR for small (below the 20th percentile), inter- 
mediate (20%— 1%) and large (above 1%) PFs over land 
and ocean, respectively. Over the oceans, the majority of 
PFs below the 20th percentile show a typical characteristic 
of shallow warm convection with storm heights at 2-4 km 
and cloud tops below or near the freezing level. Some of 
the small PFs could be as high as 8 km and colder than 
200 K and identifiable as isolated cells disconnected from 
large systems. The small PFs over the land are both higher 
and colder than their oceanic counterpart, with the highest 
frequency appearing at 5 km and 260 K. The 20%— 1 % 
range corresponds to a transition spectrum, with PFs rang- 
ing from warm and shallow systems to deep and cold 
systems. PFs in the top 1% span from middle congestus 
to deep convection, with consistent maximum height and 
minimum IR and land regions showing higher concentra- 
tions of deeper convection. Note that shallow systems 
with large VRR do occur over several specific regions of 
frequent deep convection such as the East Pacific ITCZ 
[Liu and Zipser, 2009, 2013a], 

3.7. Seasonal Variations of Extreme PF Characteristics 

[ 27 ] The large spread in the 2-D histogram (Figure 8) is 
partly due to a multitude range of precipitation systems and 
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Figure 7. Mean values with one standard deviation error bars of (a) minimum infrared brightness temper- 
ature (MinlR), (b) maximum echo top height (Maxht), (c) maximum height of 40 dBZ (Maxht40), (d) min- 
imum 85 GHz polarization correction temperature (Min85PCT), (e) minimum 37 GHz polarization 
corrected temperature (Min37PCT), and (f) flash rate for different percentiles of PFs. Note that for PFs 
without 40 dBZ echo, we used zero value in the average in Figure 7c. 


partly due to variations in the large-scale environment that 
dominates precipitation systems in different geographical 
locations and different seasons. In this subsection, we will 
examine more detailed rain characteristics of extreme PFs 
in different subregions and different seasons. For plotting 
reasons, only data from the Northern Hemisphere are used 
in the figures in this subsection (Figures 9-1 1). 

[ 28 ] Figure 9 shows the same 2-D histogram of Maxht ver- 
sus MinlR as Figure 8 but isolates the top 1% of extreme PFs 
in boreal winter (DJF) and summer (JJA) separately in four 
subregions. There are obvious differences in rain characteris- 
tics between tropics and subtropics and between winter and 
summer in the subtropics. The rain characteristics for the 
top 1% of extreme PFs in the tropics are remarkably similar 
in echo top height and cloud top temperature between the 
land and oceans (also shown in Figures 7a and 7b) and 


between summer and winter seasons, with almost complete 
overlap of the DJF and JJA contours. Extreme PFs in the 
subtropics during summer resemble those in the tropics with 
a maximum frequency of distribution located at 14 km and 
190 K, but extreme PFs in winter subtropics are significantly 
lower (6-8 km), with much warmer cloud top temperatures 
(210-220 K), which indicates very different rain systems 
such as cold fronts. 

[ 29 ] Figure 10 further shows the vertical distributions of 
raining area (represented by area with 20 dBZ radar signal) 
for the top 1% of extreme PFs. It is interesting to note 
that the largest 20 dBZ area does not occur at the surface 
but at about 2 km above the ground. The surface raining area 
is substantially smaller, with only about 30-60% of the 
area at 2 km height. There are two possible reasons for this. 
One is that the raining area identified by TRMM PR is 
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Figure 8. Two-dimensional histogram of minimum IR versus maximum stonn height for PFs in (a, b) 
<20%, (c, d) 20-1%, and (e, f) >1% categories over oceans and land, respectively. 


obtained from near-surface reflectivity measurements, 
which are valid PR observations at the lowest altitudes. 
Often, the PR profiles do not reach below 2 km because of 
ground clutter contamination. Just the lack of observations 
would lead the 20 dBZ area near the surface to be smaller, 
but a second possible reason is that evaporation of rainfall 
in the stratiform regions could also reduce the reflectivity 
near the surface. From 2 km and above, the 20 dBZ area 
decreases quickly to less than 2 x 10 4 km toward 8 km. The 
tropical rain systems extend higher vertically than the 
subtropical systems but narrower horizontally at the surface 
and near the surface. Tropical oceans have slightly lower 
vertical extension but larger rain area at lower altitudes 
than the tropical land even though maximum raining areas 
(0.1% contour) in both regions are very close in DJF. 
Tropical oceans also show an increased number of events 
with near-surface rain areas larger than 3><10 5 km 2 in 
JJA, which does not occur in the tropical land. Both 
subtropical land and oceans have small PF rainy areas at 
high altitudes (above 8 km) but have large areas at lower 
altitudes (below 5 km), with PFs in winter substantially 


larger than in summer (indicating broader and shallower 
frontal systems in winter). 

[ 30 ] The Contour Frequency Altitude Diagram (CFAD) of 
maximum dBZ shows that extreme tropical land events have 
stronger intensity, with higher maximum reflectivity values 
than those over the oceans (Figures 11a and 1 lb). At levels 
below freezing (~4.5 km in tropics), the slopes of the reflec- 
tivity increase toward the surface over the oceans but remain 
nearly vertical over land. This is possibly related to the differ- 
ence in convective intensity between land and oceans [Liu 
and Zipser, 2013b]. Tropical land has fewer contour levels 
in the center of the distribution than the tropical oceans 
(Figures 1 la and 1 lb), which is consistent with more spread 
and diversified rain system over tropical land (Figure 5). 
There is negligible difference in the vertical rain intensity 
between winter and summer in the tropics, but significant 
difference can be found in subtropical winter and summer 
(Figures 1 lc and 1 Id). Flere the vertical distributions of the 
maximum radar signals of subtropical land and oceans in 
JJA again resemble those of the tropics, but the DJF contours 
show a unique double peak feature: one with a maximum 
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four subregions. 


radar signal of 45-50 dBZ extending from near the surface to 
about 4-5 km in altitude (consistent with strong surface rain 
intensity from these extreme PFs) and the other with a max- 
imum radar signal of 25-30 dBZ located at 5-6 km above 
the ground. The latter peak should not be identified as a 
bright band because it is located above the freezing level 
(which is typically lower than 4.5 km in winter over the 
subtropics) and the maximum radar reflectivity of a PF 
usually comes from the convective region of the system 
where the bright band is not clearly detected. In fact, this 
strongly suggests dominant shallow frontal systems with a 
convective cloud top reaching only 7-9 km. 

[ 31 ] The separation of tropical and extratropical rain sys- 
tems based on geographical boundaries (20° north and 
south) will inevitably misidentify some rain systems. For 
example, some of the PFs in the tropical region may have 
moved from the extratropics and vice versa. We followed 
Xu et al. [2007], using an eddy growth parameter o as a 
baroclinic instability criterion to physically separate the 
tropical and extratropical rain systems. The eddy growth 
parameter is computed as o = 0.31 x f/N x (d|V|/dZ), where 
f is the Coriolis parameter, N the Brunt-Vaisala frequency, 
and djV|/dZ is computed from winds at 1000 and 700 hPa. 
If o is greater than 0.1 /day, the system is considered to 
be baroclinic and extratropical in nature; otherwise, it is 


considered tropical. As a result, we found that about 10% 
and 18% of PFs in DJF in subtropical land and oceans, re- 
spectively, may belong to tropical systems because of their 
small baroclinic instability. In JJA, nearly half of the PFs 
do not meet the baroclinic instability criterion. However, 
the structural difference with respect to the baroclinic crite- 
rion is rather insignificant in the subtropical regions. In the 
tropical regions (0°N-20°N), about 20-30% of PFs meet 
the baroclinic criterion and may be considered to be 
extratropical systems. We notice that the PFs with larger 
baroclinic instability in the tropics are relatively larger in 
size and have lower maximum dBZ levels (figures not 
shown), which are consistent with the overall results from 
Figures 10 and 11. 

[ 32 ] The above results (Figures 9-11) show that extreme 
PFs in the tropics extend higher with stronger rain intensity 
than the subtropical rain systems and experience very little 
seasonal variation. Summertime extreme PFs in the subtrop- 
ics, characterized by smaller rain areas but thicker and 
heavier vertical rain volumes, resemble those in the tropics 
and is consistent with more locally convective-driven sys- 
tems. Wintertime subtropical extreme PFs are lower and 
broader, indicating characteristics of frontal rain systems. 
Of all the subregions studied, seasonal variation is most 
significant in subtropical land. 
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Figure 10. Cumulative frequency with altitude (CFADs) of the area with 20 dBZ for the top 1% of ex- 
treme events in DJF (shading) and JJA (contour) in four subregions. 


4. Large-Scale Environment 

[ 33 ] Many studies have attempted to characterize typical 
meteorological conditions and major patterns of precipitation 
systems associated with extreme precipitations on a regional 
basis [. Maddox et al., 1979; Glass et al., 1995; Doswell et al., 
1996; Moore et al., 2003; Schumacher and Johnson, 2005; 
Romatschke and Houze, 2010, 2011a, 2011b]. The synoptic 
and mesoscale environments of extreme precipitation from dif- 
ferent geographical regions and for different storm types can be 
quite varied. Here we composite the large-scale environments 
for different percentiles of PFs to see if such systematic 
changes exist. We show 500 hPa relative humidity, total pre- 
cipitable water, 500 hPa vertical velocity, vertical wind shear, 
and convective available potential energy (CAPE) in particular 
to represent atmospheric thennodynamic, dynamic, and con- 
vective potentials, respectively. 

[ 34 ] Figure 12a shows that mean precipitable water re- 
mains constant for PFs below the top 10th percentile but 
increases rapidly after the PFs reach the 10th percentile and 
beyond, suggesting that increased total precipitable water is 
essential to extreme rainfall. There is a clear separation of 
the total precipitable water between the tropics and subtrop- 
ics regardless of land or ocean, with the tropics being about 
1 5 mm higher than the subtropics due to wanner temperature 
and higher specific humidity. The differences narrow as the 


mean precipitable water in all regions reaches around 50 mm 
for the extreme PFs. It is interesting to notice that the 
500 hPa mean relative humidity over land is slightly higher 
than the relative humidity over the oceans for small PFs 
between the 90th and 1 0th percentiles, indicating that for small 
precipitations to occur, the moistening of the free atmosphere 
is even more important over land than over the oceans. The 
mean 500 hPa relative humidity rises rapidly with rainfall in 
all the regions until the PF reaches the critical threshold of 
the 10th percentile but levels off for the most extreme ones 
in the subtropics (Figure 12b). The surge of midlevel relative 
humidity for extreme PFs is likely related to the large-scale 
moisture convergence and convective uprising that brings in 
additional moisture to that level. 

[ 35 ] The rapid increase in relative humidity for the extreme 
PFs is well correlated with an increase in 500 hPa vertical 
velocities. Similar to relative humidity, the mean 500 hPa 
vertical velocity remains constant at around 0. 1 Pa/s for small 
PF s up to the 1 0th percentile, then increases abruptly after the 
10th percentile and reaches tenfold in magnitude in the top 
0.01% compared to those in the lower 90% (Figure 12c). 
For extreme PFs, subtropical storms have higher vertical 
velocity compared to tropical stonns, which might be due 
to stronger convergence and vertical motions, both required 
to support the required moisture convergence in the relatively 
dry area (less total precipitable water in subtropics). It is well 
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Figure 11 . Contoured frequency by altitude diagrams (CFADs) of maximum PR reflectivity for the top 
1% of extreme events in DJF (shading) and JJA (contour) in four subregions. 


known that precipitation at any given location comes from 
evaporation and moisture convergence. For heavy precipita- 
tion, moisture convergence is far more important than local 
evaporation. The rapid increase of midlevel relative humidity 
and vertical velocity shows that large-scale dynamics plays a 
critical role to local precipitation [O 'Gorman and Schneider, 
2009; Chou et al., 2012]. 

[36] Besides atmospheric moisture and vertical velocity, 
which are directly related to extreme precipitation, many 
other environmental factors could affect the generation, 
intensity, and duration of precipitation systems. For example, 
large environmental vertical wind shear hinders tropical 
cyclone development but favors multicell systems in the 
subtropics. Strong ambient vertical shear is found to be 
very important to the structure, strength, and longevity of 
subtropical squall lines over land [Bluestein and Jain, 
1985; Rotunno et al., 1988; Weisman and Rotunno, 2004], 
To examine whether vertical wind shear affects rain volume 
in any systematic manner, we plotted mean vertical shears 
between 250 and 850 hPa for different percentiles of PFs 
(Figure 12d). It is not surprising to see that mean vertical 
wind shear values over the subtropics are much higher than 
those over the tropics due to a more baroclinic atmosphere 
in general. Vertical shear varies little for the small PFs for 
all the regions. Flowever, they start to increase at the 10th 
percentile in subtropical land and oceans, indicating larger 
vertical shear is conducive to some subtropical stonn systems. 


Vertical wind shear in both tropical land and oceans remains 
constant for the entire rain spectrum even with the known 
effects on tropical cyclones, probably due to the small percent- 
age of tropical cyclone-related PFs of the total PFs. 

[ 37 ] The last environmental parameter examined is related 
to the atmospheric instability, i.e., convective available 
potential energy (CAPE). Since convective related precipita- 
tion is due to buildup of atmosphere instability, precipitation 
through vertical overturning and latent heat release will 
remove or reduce that instability. Figure 12e shows large dif- 
ferences in the mean CAPE from different regions; the mean 
CAPEs over subtropical land and oceans are much smaller 
than those of the tropics due to a large number of subtropical 
systems having very small or no CAPE. The CAPE generally 
increases slightly with PF percentiles up to 10% for all the 
regions but only continues increasing slowly for extreme 
PFs in tropical oceans. For the subtropical land region, there 
is a discernible decreasing trend of CAPE from around 850 J/ 
kg at the top 10% to around 650 J/kg at the top 0.01%, which 
is likely due to the release of some of the prestorm instability 
in these large mature systems. Flowever, it should be noted 
that CAPE is a loose measure of atmospheric instability and 
the standard deviation is very large (~700 J/kg). How effec- 
tive the instability can be converted into actual convective 
activity (i.e., vertical velocity and rainfall) depends on many 
factors such as large-scale moisture convergence, convective 
inhibition, midlevel moisture, wind shear, and entrainment 
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Figure 12. Mean values with one standard deviation error bars of (a) 500 hPa relative humidity, (b) total 
precipitable water vapor, (c) 500 hPa vertical velocity, (d) 200-850 hPa wind shear, and (e) CAPE from 
NCEP reanalysis for PFs with different percentile of volumetric rainfall. 


[Lucas et al., 1994]. Therefore, we should expect a generally 
weak relationship between VRR and CAPE and that the same 
CAPE would work differently over different regions. In addi- 
tion, due to the temporal interpolation and coarse horizontal 
resolution, there is a large uncertainty in the CAPE computed 
from the NCEP reanalysis, especially for small PFs. As a 
result, these relationships should be considered qualitative 
in nature. 

[38] This section shows very nonlinear behavior of the 
large-scale parameters that experience little change for the 
lower 90% of PFs but change rapidly after the 10th percen- 
tile. These environmental conditions are directly related to 
the observed abrupt changes in the structure of the PF, which 
is likely caused by convective-driven vertical motion accom- 
panied by large-scale moisture convergence and the release 
of CAPE. 

5. Summary and Discussions 

[ 39 ] This study adopts a “precipitation object” approach by 
using the Precipitation Feature (PF) database developed at 
the University of Utah to analyze the characteristics and 
large-scale environment of extreme precipitation from 14 year 
TRMM observations. We define extreme PF with the volu- 
metric rain of the PF in order to capture the effect of both 


spatial extent and rain intensity. We have found that instanta- 
neous PFs have much larger dynamic ranges than the daily 
gridded precipitation, with those in the top 1% being 2 orders 
of magnitude larger than the medium PFs and contributing to 
over 55% of instantaneous rainfall. Ninety percent of the PFs 
that contribute about 20% of total precipitation belong to 
meso-y systems (less than 20 km in the horizontal direction) 
or broken patches from large systems. The most abrupt 
change is found to be around the top 1 0% — 1 %, where the rain 
systems expand from a few hundred to a few thousand kilo- 
meters in size and from 2300 to 40,000 mm/h km 2 in rain 
volume. Spatial and temporal distributions and the mean 
characteristics of extreme PFs are discussed in comparison 
to intermediate and smaller PFs. We find that regional differ- 
ences in rain characteristics mainly come from land-ocean 
differences for smaller PFs. As the VRR of PF increases, 
the regional differences mainly come from the tropics and 
subtropics. While extreme PFs in the top 1% are found to 
be significantly larger, deeper, and colder than the lower 
90% of PFs, extreme PFs in the tropics are even deeper and 
colder than the subtropical systems and present no significant 
seasonal variations. The extreme PFs in the subtropics in the 
summertime resemble those in the tropics, characterized by 
locally convection-driven rain systems. The extreme PFs in 
subtropical winter are much broader and shallower, which 
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is a signature of midlatitude frontal systems. The NCEP 
reanalysis shows that midlevel relative humidity and total 
precipitable water increase consistently with increasingly 
larger rain systems, along with sharp increases of upward 
vertical velocities to provide necessary dynamic support. 
Instantaneous CAPE also generally increases with larger 
PF, but the mean CAPE levels off for the extremely large 
systems due to the release of CAPE in mature systems in 
subtropical land regions. The study illustrates important roles 
of large-scale moisture and dynamic conditions for occur- 
rences of extreme precipitation. 

[ 40 ] It has been widely accepted that increasing moisture 
in the atmosphere potentially provides a mechanism for 
increasing extreme precipitation in a warming climate 
[Trenberth et al., 2003; Allan and Soden, 2008]. However, 
the correlation between surface temperature, atmospheric 
water vapor, and precipitation is not always positive and is 
highly dependent on the regions and time scales [Adler 
et al., 2008; Gu and Adler, 201 1], In addition, both observa- 
tions and climate model simulations show that increases in 
extreme precipitation do not always follow the rate of 
increase in atmosphere water content [Lau and Wu, 2011]. 
Therefore, it is important to confirm that the total water 
vapor in the atmosphere is indeed highly correlated to 
the total volumetric rainfall, even on an instantaneous 
basis. This proves that extreme rainfall is more likely to 
occur when high water vapor content is readily available 
in the atmosphere. 

[ 41 ] This study, as well as previous studies [Doswell 
et al., 1996; Schumacher and Johnson, 2005], points out 
the necessary conditions for extreme precipitation to occur, 
such as high atmospheric water vapor, high upward vertical 
velocity, and high atmospheric moist static energy. It is 
tempting to consider that each of these meteorological 
parameters can reach their optimum statuses when a preci- 
pitation system is organized, which might have no direct 
connection to global warming. However, it has been shown 
that it is the thermodynamic effect rather than the dynamic 
effect that is considered mainly responsible for changes 
in precipitation intensity [Allan and Soden, 2008; Emori 
and Brown, 2005; O’Gorman and Schneider, 2009; Chou 
et al., 2012], Thus, water vapor in the atmosphere is the 
most relevant parameter related to the effect of global 
warming. This study shows a consistent increase of mean 
volumetric rainfall with the increase of total precipitable 
water in the real world, which could imply that there is more 
extreme rain in a wanning climate when atmospheric water 
vapor loading is high. 

[ 42 ] Obviously, there are large standard deviations for the 
mean characteristics and large-scale environments for all 
percentiles of the PFs. This is not a surprise because many 
different types of rain systems occur in each location and 
nature never repeats itself. Therefore, these results should 
be interpreted as qualitative rather than quantitative. The 
other caveat is related to the instantaneous nature of PFs; 
some of the extreme PFs do not produce extreme rainfall 
accumulation. However, our results are based on a large 
number of PF samples over the global tropics; thus, the 
observed systematic changes of the rainfall characteristics 
and large-scale parameters associated with extreme PFs are 
fundamental to understanding the mechanisms of develop- 
ment of extremes in current and future climates. 
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